from Excitation_Carbon import *

cs = CS_Excitation(species_name = "C+1")


id = 0

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 4.74,
    A   = 4.360e-1,
    B   = 1.627,
    C   = -2.923,
    D   = 1.094e1,
    E   = -6.580,
    F   = 5.769e-2,
    omega = 6.0,
    equation_number = 7,
    reaction_name = "Excitation_C+1_2s22p2P-2s2p24P")


id = id + 1
cs.add_cross_section(
    id = id,
    Vif = 9.27,
    A   = -2.155,
    B   = 1.374e1,
    C   = -6.562,
    D   = 0.0,
    E   = 1.018,
    F   = 0.0,
    omega = 6.0,
    equation_number = 6,
    reaction_name = "Excitation_C+1_2s22p2P-2s2p22D")


id = id + 1
cs.add_cross_section(
    id = id,
    Vif = 12.0,
    A   = 8.572e-1,
    B   = 2.079e-1,
    C   = 7.058e-1,
    D   = 0.0,
    E   = 4.136,
    F   = 0.0,
    omega = 6.0,
    equation_number = 6,
    reaction_name = "Excitation_C+1_2s22p2P-2s2p22S")

        


id = id + 1
cs.add_cross_section(
    id = id,
    Vif = 13.7,
    A   = -6.883e-1,
    B   = 4.106,
    C   = 3.435e-1,
    D   = 0.0,
    E   = 1.631e1,
    F   = 0.0,
    omega = 6.0,
    equation_number = 6,
    reaction_name = "Excitation_C+1_2s22p2P-2s2p22P")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 14.4,
    A   = -9.843e-1,
    B   = 2.983,
    C   = -9.975e-1,
    D   = 0.0,
    E   = 8.910,
    F   = 0.0,
    omega = 6.0,
    equation_number = 6,
    reaction_name = "Excitation_C+1_2s22p2P-2s23s2S")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 16.3,
    A   = 1.678,
    B   = -2.238,
    C   = 4.829,
    D   = 0.0,
    E   = 0.0,
    F   = 0.0, 
    omega = 6.0,
    equation_number = 6,
    reaction_name = "Excitation_C+1_2s22p2P-2s23p2P")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 18.0,
    A   = -1.703,
    B   = 2.675,
    C   = 1.165,
    D   = 0.0,
    E   = 5.791,
    F   = 0.0,
    omega = 6.0,
    equation_number = 6,
    reaction_name = "Excitation_C+1_2s22p2P-2s23d2D")

cs.plot_cs()
cs.save_cs()
cs.cal_total_cs()